The development of generalized synchronization on complex networks 
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In this paper, we investigate the development of generalized synchronization (GS) on typical complex net- 
works, such as scale-free networks, small-world networks, random networks and modular networks. By adopt- 
ing the auxiliary- system approach to networks, we show that GS can take place in oscillator networks with 
both heterogeneous and homogeneous degree distribution, regardless of whether the coupled chaotic oscillators 
are identical or nonidentical. For coupled identical oscillators on networks, we find that there exists a general 
bifurcation path from initial non- synchronization to final global complete synchronization (CS) via GS as the 
coupling strength is increased. For coupled nonidentical oscillators on networks, we further reveal how net- 
work topology competes with the local dynamics to dominate the development of GS on networks. Especially, 
we analyze how different coupling strategies affect the development of GS on complex networks. Our findings 
provide a further understanding for the occurrence and development of collective behavior in complex networks. 

PACS numbers: 05.45.Xt,89.75.Hc 



I. INTRODUCTION 

Synchronization in coupled chaotic oscillators has been 
extensively studied in the past twenty years [1], including 
complete synchronization yfl, generalized synchronization 
HHSBEI], and phase synchronization (PS) |8], etc. In CS, 
the dynamics of two coupled systems totally coincide with 
each other; in GS, certain functional relation exists between 
the dynamics of two coupled systems which are usually non- 
identical. Moreover, PS is a weaker synchronization form 
in which the phases of two oscillators can be locked while 
their amplitudes remain uncorrelated and chaotic. Recently, 
the study of synchronization has been extended to the area of 
complex networks M O d H H H M [HI MM 
I^l2lil22ll231. F or example, synchronization on small-world 
networks l9l 1 1 0L fill fl2|j. scale-free networks HH Q3, fl5[l . 
modular networks 1 16, 17, , weighted networks [ 20], 

and gradient networks 12111 . have been investigated. These 
studies aim to explore the interplay between network topol- 
ogy and dynamics on network. They are important for us to 
understand the real situations in complex systems comprising 
interacting elements in nature and human society. 

In literatures, there are two main approaches studying syn- 
chronization in complex networks. One approach is to study 
synchronization of phase oscillators on complex networks, 
i.e., the generalized Kuramoto models [S [II S S Hi- 
In this model, the node dynamics is very simple which is 
governed by an ideal phase oscillators, (j) = u, where uj is 
the frequency. The heterogeneity in the node dynamics can 
be modeled by assigning different, usually random, frequen- 
cies to different phase oscillators. The generalized Kuramoto 
models in complex networks have the advantage that they can 
still be treated analytically in many aspects The sec- 

ond approach is to study the synchronizability of complex net- 
works through the master stability function j 111 12511 . To apply 
this approach, the node dynamics in complex networks are as- 
sumed to be identical, and then the theory of master stability 
function provides a general mathematical framework to relate 



the synchronizability of a network to the spectral properties of 
the corresponding coupling matrix lllll[l4flAl^ l26n. 

The above two approaches allows an analytical and there- 
fore relatively more complete understanding of synchroniza- 
tion in complex networks. Nevertheless, in realistically physi- 
cal or biological situations the node dynamics can conceivably 
be much more complicated, and in general be different from 
each other. The key point is that the possible collective behav- 
iors for coupled nonidentical oscillators should be of the type 
of GS rather than CS and PS. It is thus of great importance 
to investigate the collective behaviors on complex networks 
where the local node dynamics might be different from each 
other. In this paper, we present a systematical study on GS 
phenomenon in various complex networks. Our particular at- 
tention is paid to the development of GS in networks. For 
typical complex networks, such as scale-free networks, small- 
world networks, random networks and modular networks, in- 
terestingly, we find that GS generally occurs in these net- 
works, regardless of whether the node dynamics are identical 
or nonidentical. For coupled identical oscillators on complex 
networks, a general bifurcation path toward the final global 
CS is shown as: non- synchronization — » PGS — » GS — » PCS 
-> CS. Here, PGS and PCS stand for partial GS and partial 
CS, respectively. We further study how the network topology 
competes with the local dynamics to affect the development 
of GS on networks. It is shown that for linearly coupled oscil- 
lators on heterogeneous networks, GS usually starts from the 
hub nodes and then spreads to other nodes; while on homo- 
geneous networks, GS generally starts from the nodes whose 
local dynamics are less chaotic in the sense that the largest 
Lyapunov exponents have relatively smaller values. More- 
over, we analyze how different coupling strategies can affect 
the development of GS on complex networks, a point which 
has not been addressed before. 

In processing this paper, we noticed a recent work which 
has reported GS phenomenon in scale-free networks Il27n . 
Compared with the results in Ref. B27M . the present work is 
different in the following aspects. Firstly, in our work GS is 
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systematically investigated on various complex networks, in- 
cluding scale-free networks, random networks, small- world 
networks and modular-networks; while in Ref. Il27ll . GS is 
mainly studied on a very special network, i.e., the scale- 
free networks with tree-like structure. Secondly, the present 
work characterizes the complete bifurcation path for cou- 
pled identical oscillators on complex networks, i.e., from 
non- synchronization to global CS, including the transition 
regimes; while in Ref. 112711 . global CS has not be achieved. 
Thirdly, the present work investigates the development of GS 
for coupled nonidentical oscillators, either parametrically dif- 
ferent or physically different; while Ref. [21] does not con- 
sider this general setting where the occurrence of GS is nat- 
urally justified. Finally, in the present work, we discuss the 
effect induced by different coupling strategies on the develop- 
ment of GS on network. Especially, we analyze and explain 
the different observations regarding the GS evolution on net- 
works in this work and in Ref. Il27ll . In fact, in many aspects 
the current study deepens and widens the work in Ref. ll27h . 
and can offer more thorough and comprehensive insight in un- 
derstanding GS phenomenon on complex networks. 

This paper is organized as follows. In Sec. II, the meth- 
ods to characterize GS and CS on complex networks are de- 
scribed. In Sec. Ill, GS of coupled identical oscillators is 
studied on complex networks. The complete bifurcation path 
from non- synchronization to global CS is studied and the dy- 
namical features in the transition regimes are characterized. 
In Sec. IV, GS of coupled nonidentical oscillators is further 
studied on complex networks. Our interest is focused on how 
the development of GS on networks can be affected by net- 
work topology and the local dynamics. In Sec. V, the effect 
of different coupling strategies on the development of GS on 
complex networks is analyzed. Especially, GS is observed in 
coupled system with hybrid oscillators even when their local 
dynamics are physically different. Finally, a summary ends 
this paper in Sec. VI. 

II. APPROACHES CHARACTERIZING GS AND CS ON 
NETWORKS 

The auxiliary-system approach has been extensively used to 
detect GS in two coupled chaotic systems |5]. Here, we can 
extend it to detect GS on complex networks. The key obser- 
vation is that, for any given node in a network, the coupling 
from other nodes can be regarded as a kind of "driving." In 
particular, we consider the following linearly coupled identi- 
cal oscillators on a network: 

±i Fifa) -s^aij^i - Xj), (1) 

3 

for i = 1, . . . , N, where denotes the dynamical variables 
of node i, F^x^) is the local vector field governing the evolu- 
tion of x^ in the absence of interactions with other nodes, 
is the element of the network adjacency matrix A (a^ = 1 
if there is a link between node i and node j, = other- 
wise, and an = 0), and e is the coupling strength. To apply 
the auxiliary-system approach, we consider a replica for each 



oscillator in the original network: 

x^F^xO-e^a^-x,-), (2) 

3 

for i = 1, . . . , N. Note that the driving variable x^ is iden- 
tical for both Eqs. (Q} and 0. If, for initial conditions 
x^(0) 7^ x-(0), we have |x^(t) — x-(t)| — > as t — > oo, node 
i then is entrained in the sense that its dynamics is no longer 
sensitive to initial conditions. In other words, there is GS re- 
lation between x^ and Xj for j — 1, . . . , N. Numerically, 
we can examine the following local distance of GS between a 
node and its auxiliary counterpart: 

1 t2 

d(e,i) = - r £ l*i I > (3) 

ti 

where ti is chosen to be larger than the typical transient time 
of the local dynamics F^x^). For oscillators on complex net- 
works, GS may be gradually developed with the increase of 
the coupling strength. To characterize the development of 
GS on networks, we can define the distance of global GS as 
l g (e) = (d(e,i)). Here, (•) denotes the spatial average over 
all nodes. If l g = 0, global GS has been achieved between any 
two pair of oscillators on the whole network. 

For coupled identical oscillators on complex networks, CS 
is generally expected to take place as long as the network con- 
nectivity is not too sparse. To characterize CS, we can define 
l c (e) as the distance of global CS, which measures the dis- 
tance between the dynamics of all oscillators and their aver- 
age, i.e., 

^) = ^^X>-<x>i>, (4) 

ti 

where the meaning of t\ and is the same as that in Eq. ©. 
If l c = 0, global CS has been achieved on the whole network. 

For individual node dynamics, we choose the chaotic 
Lorenz oscillator: 

Fi(xi) = [10(xi-yi),riXi-yi-XiZi,Xiyi-(8/3)zi] T , (5) 

where x^ = (a^, y%,Zi) are the state variables of the Lorenz 
oscillator. Note that in studying GS, the local dynamics could 
be different from each other. This can be modelled by setting 
different T\ values for Lorenz oscillators in the network. With- 
out losing generality, the coupling between two nodes in Eq. 
(Q} is through the x variable. For convenience, in the present 
work, we use two kinds of node indices: id and i r . In the first 
index, we order the degrees of the network so that id = 1 de- 
notes the node with the largest degree, id — 2 is for the node 
with the second largest degree, and so on. In the second in- 
dex, we rank the Lorenz oscillators in the network according 
to their n values, i.e., i r = 1 denotes the node with the largest 
parameter r, i r = 2 is for the node with the second largest 
parameter r, and so on. These two indices of nodes are used 
throughout the paper. 
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FIG. 1: (Color online) Color-map of d(s, i) in the two-dimensional 
parameter space (e, i), characterizing the development of GS for 300 
identical chaotic Lorenz oscillators on a scale-free network. The net- 
work is the BA model with mo — m — 4 [13]. We see that with the 
increase of coupling strength, GS can be achieved gradually from the 
hubs and then spread to the rest of the network. 



III. COUPLED IDENTICAL OSCILLATORS: A GENERAL 
PATH TOWARD GLOBAL CS 

A. Path towards CS on scale-free networks 
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FIG. 2: (Color online) Characterizing the path towards global CS on 
scale-free network, (a) The distance of global GS vs the coupling 
strength, (b) The distance of global CS vs the coupling strength. 
Numerically, if l g (l c ) < 0.001, global GS (CS) is regarded to have 
been achieved. The notations of the regimes in this figure are used 
throughout this paper. 



Previously, in the study of synchronization of two coupled 
chaotic oscillators, it is found that CS generally takes place 
when the coupled oscillators are identical, while GS usually 
occurs in systems of coupled nonidentical oscillators. This 
observation is accurate for most cases except for some spe- 
cially designed coupling strategies |2D- However, this gen- 
erally accepted view turns out to be not true when synchro- 
nization is studied on complex networks. In this section, we 
investigate the occurrence of GS for coupled identical chaotic 



oscillators on various complex networks. Our particular inter- 
est is focused on the complete bifurcations from initial non- 
synchronization towards the final global CS with the increase 
of the coupling strength. By applying auxiliary-system ap- 
proach, interestingly, we find that GS generally exists in cou- 
pled identical oscillators on various complex networks, in- 
cluding scale-free networks, random networks, small- world 
networks and modular networks. 



As an example, we first study the bifurcation route towards 
global CS for 300 coupled chaotic Lorenz oscillators on a 
scale-free network. In this case, r^=28 in Eq. ([5]) for all nodes, 
i.e., all oscillators are identical. Intuitively, with the increase 
of coupling, we expect that the coupled chaotic oscillators will 
finally achieve global CS. The surprising finding here is that 
before the system achieves the final global CS, there exists 
another synchronization regime, i.e., the GS regime, which 
usually occurs with much smaller coupling strength compared 
with CS. To illustrate the GS phenomenon on scale-free net- 
work, we plot the color-map of the distance matrix d(s,id) 
in Fig. 1. From Fig. 1, three distinct regimes of coupling 
can be identified from the viewpoint of GS. When the cou- 
pling strength is small, d(e, id) are greater than for all nodes, 
showing that the system is in the non- synchronous state. On 
the other hand, when the coupling strength is large enough, 
d(e, id) are for all nodes, showing that all oscillators have 
been entrained and the coupled system is in the global GS 
state according to the auxiliary system approach criterion. Be- 
tween these two regimes, it is the regime of partial GS, where 
part of the oscillators have been entrained but the others have 
not. Obviously, this regime is a transition stage between non- 
synchronous state and the global GS state. 



To quantitatively characterize the GS regime for coupled 
chaotic oscillators on scale-free network, in Fig. 2(a) we plot 
the distance of global GS l g vs the coupling strength. From 
this figure, three regimes of coupling can numerically be iden- 
tified. On the other hand, for coupled identical oscillators, CS 
is generally expected to occur as long as the connections of 
the network are dense enough. It is found that usually CS re- 
quires larger coupling strength than GS. This implies that GS 
is a weaker synchronization form compared with CS on com- 
plex networks, as in the situation of coupled two-oscillator 
systems. To characterize CS in the network, we plot the dis- 
tance of global CS l c vs the coupling strength in Fig. 2(b). 
We find that before CS is achieved, there is also a transition 
stage, i.e., the partial CS regime. Combining all the results in 
Fig. 2, we can identify five dynamical regimes on the path to- 
wards global CS on the specific scale-free network as follows. 
lie < 0.1, the non- synchronization regime; II: 0.1 < £ < 1.8, 
the PGS regime; III: 1.8 < e < 3.0, the global GS regime; IV: 
3.0 < e < 3.8, the PCS regime; V: £ > 3.8, the global CS 
regime. These regimes thus characterize the complete bifur- 
cations from non- synchronization state to global CS state for 
coupled identical oscillators on scale-free networks. 
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FIG. 3: (Color online) Typical routes toward final GS for Lorenz 
oscillators in scale-free network. From (a) to (c), d(e,id) versus e 
for id =4, 104, and 250, respectively. We observe that nodes usually 
oscillate between the entrained state (GS) and the unentrained state 
on the way toward final global GS state as the coupling strength is 
increased. 





FIG. 5: (Color online) A schematic visualization of the oscillatory 
path toward global GS shown in Fig. 4(a). (a) Non- synchronous 
state at e = 0. (b) PGS state at si = 0.55. About 80 percent nodes 
in the network (nodes with blue color) are entrained, (b) PGS state at 
si — 0.61. About 50 percent nodes in the network (nodes with red 
color) are entrained, (d) Global GS state at e = 3. 
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FIG. 4: (Color online) Characterizing the dynamic and oscillatory 
nature in the PGS and PCS regimes on a scale-free network, (a) The 
number of nodes which have achieved GS vs the coupling strength, 
(b) The number of nodes which have achieved CS with the hub (the 
node with the largest degree) vs the coupling strength. 



B. Dynamic and oscillatory transition to global 
synchronization 

In the above we have shown that there are two transition 
regimes, i.e., the PGS regime and PCS regime, before the cou- 
pled system achieves global GS and CS, respectively. In fact, 
these two regimes are important to obtain information about 
how synchronization is generated and spread on a network. 
As shown in Fig. 1, with the increase of coupling strength, 



usually the hubs of scale-free network become entrained first. 
Then more and more nodes are entrained with further increase 
of the coupling strength. At £ « 0.3, a giant cluster of syn- 
chronous nodes in the sense of GS has been formed. The most 
important dynamic feature in the PGS regime is that nodes 
may oscillate between the entrained state and the unentrained 
state with the increase of the coupling strength. As shown 
in Fig. 3, if we focus on a specific node, we find that its 
path toward global GS is generally oscillatory. Especially, for 
those nodes with small degrees, the path could be extremely 
oscillatory as shown in Fig. 3(c). As a consequence, in the 
PGS regime, the size of the synchronous giant cluster does 
not increase monotonically, instead it increases in an oscil- 
latory manner with the increase of the coupling strength, as 
shown in Fig. 4(a). Similar phenomenon has been observed 
in the PCS regime. For example, in Fig. 4(b), we plot the 
number of nodes which have achieved CS with the hub (with 
the largest degree). This figure actually characterizes the dy- 
namical features for regime IV shown in Fig. 2(b). It is found 
that the size of the synchronous cluster oscillates drastically 
on the way toward global CS with the increase of the coupling 
strength. 

Recently, the path towards global PS in generalized Ku- 
ramoto model has been studied in Ref. 11231] . It is found that 
for phase oscillators with the simplest local dynamics, once a 
phase oscillator becomes frequency locked, it will never es- 
cape from the synchronous cluster. Thus in the PS situation 
the size of the synchronous cluster increases monotonically 
with the increase of the coupling strength. Such a path toward 
global synchronization can be called as static and monotonic. 
In the present work, we demonstrate a fundamentally different 
scenario towards global synchronization for coupled chaotic 
oscillators on complex networks. The most important feature 
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of the path is that the development and spreading of synchro- 
nization on networks is a dynamic process during which oscil- 
lators may frequently oscillate between the non-synchronous 
state and the synchronous state. This implies that the size of 
the synchronous cluster might be sensitive to the perturbation 
of the coupling strength in the transition regimes. To visualize 
the dynamic feature of the development of synchronization on 
networks, a schematic figure has been plotted in Fig. 5 for a 
scale-free network. Our results suggest that in more realistic 
situations, the path towards global synchronization on com- 
plex networks could be dynamic and oscillatory rather than 
static and monotonic. 



100 

3 

200 
300 

i 

300 
D 200 
100 


I 

300 
200 

i 

100 





I 


WW m 

f • 


in 


0.5 1 


1.5 2 






ill lv 


V 



(b) 



(c) 



2 3 

e 



FIG. 6: (Color online) Characterizing the development of GS and the 
path toward global CS of 300 coupled identical Lorenz oscillators 
on a random network with average degree 10. (a) The color-map 
of d(£,i) in the two-dimensional parameter space (b) The 

number of nodes which have achieved GS vs the coupling strength, 
(c) The distance of global CS vs the coupling strength. 



case of complex networks, a regular network consisting of 100 
nodes and each node has 10 nearest neighbor connections. 

The pathes towards global CS on these typical complex 
networks have been characterized as for scale-free networks. 
Here, we only show one example in the case of random net- 
work. In Fig. 6, the bifurcation route to global CS for 300 
identical Lorenz oscillators on a random network is character- 
ized. Comparing Fig. 6(a) with Fig. 1, we can see that the de- 
velopment of GS on random network is different from the situ- 
ation on scale-free network. Mainly, most nodes are entrained 
at approximately the same coupling strength. This is consis- 
tent with the fact that random network has approximately ho- 
mogeneous degrees and thus no extraordinary hubs exist as 
in scale-free networks. Although the development of GS on 
different networks depends on the network topologies, qual- 
itatively, the path toward global GS is the same as shown in 
Fig. 1 and Fig. 6. In fact, for other types of complex networks, 
we have similar observations, as in the scale-free network and 
random network. Based on these studies, we can draw the fol- 
lowing conclusions. (1) For coupled identical chaotic oscilla- 
tors, whether the network is homogeneous or heterogeneous, a 
general bifurcation path towards global CS exists with the in- 
crease of coupling strength, i.e, non- synchronization — > PGS 
-> GS -> PCS -> CS. (2) Usually, global GS needs relatively 
smaller coupling strength to achieve compared with global 
CS. Sometimes, the latter requires ten times larger coupling 
strength to achieve than the former. This shows that GS is a 
weaker synchronization form on complex networks compared 
with CS, and is consistent with the observation in coupled 
two-oscillator systems. (3) Typically, the coupled system un- 
dergoes a dynamic and oscillatory transition to achieve global 
synchronization, such as global GS or global CS. During the 
PGS or PCS regime, oscillators may frequently join in and es- 
cape from the giant synchronous cluster. This leads to that the 
size of synchronous giant cluster increases in an oscillatory 
manner with the increase of the coupling strength. In some 
circumstances, the transition regimes cover quite wide range 
of coupling strength, reflecting that there exists intense inter- 
action among nodes in achieving global synchronization on 
complex networks. 



C. Path toward global CS on other types of complex networks 

In the above we have characterized the path towards global 
CS for coupled identical chaotic oscillators on scale-free net- 
work, which is characterized by a heterogeneous degree dis- 
tribution. We have also investigated synchronization path for 
coupled identical chaotic oscillators on other typical complex 
networks. Specifically, the networks studied include the fol- 
lowing: (1) A random network consisting of 300 nodes. The 
average degree is 10. (2) A small- world network, which is ob- 
tained by rewinding 20 links in a regular network consisting of 
100 nodes and each node has 10 nearest neighbor connections 
|@]. (3) A modular network consisting of 100 nodes which 
is evenly divided into 5 modules. Inside each module nodes 
are fully connected. Any two nodes in different modules have 
probability p = 0.01 to connect each other. (4) As a special 



IV. COUPLED PARAMETRICALLY DIFFERENT 
OSCILLATORS: TOPOLOGY VS LOCAL DYNAMICS 

So far, most of the works on synchronization in complex 
networks study the coupled identical oscillators, or coupled 
simple phase oscillators of Kuramoto type. However, in real- 
world networks, the interacting components are in principle 
different from each other. For example, in neuron networks 
any two pair of neurons cannot be exactly the same. This 
raises an important question: is it possible for complex net- 
works with different local dynamics achieve GS, or certain 
coherence? If so, how network topologies affect the GS? In 
this section, we will investigate synchronization of coupled 
nonidentical oscillators on complex networks. Our particular 
interest is to reveal how the development of GS on networks 
is dominated by the interplay between network topology and 
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FIG. 7: (Color online) Color-map of d(e, id) in the two-dimensional 
parameter space (e, i), characterizing the development of GS for 500 
non-identical Lorenz chaotic oscillators on a scale-free network. The 
network is the BA model with mo = m = 3 fl3ll . The development 
of GS in the network is similar to the situation shown in Fig. 1. 



the local dynamics. 

As in the case of coupled system with two different oscilla- 
tors, we find that CS is generally forbidden for coupled non- 
identical oscillators on complex networks. Instead, we ob- 
serve that GS generally occurs in such coupled systems on 
networks. In our study, the Lorenz oscillators, i.e., Eq. (0), 
are made different on network by randomly setting parameter 
Ti in the interval [28.0,30.0]. By applying auxiliary system 
approach, we can plot the color-map of the local synchroniza- 
tion distance d(e, id) in the two-dimensional parameter plane 
(e,id). As an example, Fig. 7 illustrates the development 
of GS for Lorenz oscillators with parameter mismatches on a 
scale-free network. It is shown that with the increase of cou- 
pling strength, nodes in the network are entrained first from 
the hubs and gradually spread to other nodes. When the cou- 
pling strength is large enough, all nodes in the network are 
entrained indicating global GS has been achieved. This sit- 
uation is similar to the case of coupled identical oscillators 
due to the fact that the parameter mismatches here are small. 
Moreover, for coupled nonidentical oscillators, the dynamical 
feature of the route toward global GS is essentially the same 
as in the case of identical oscillators. 

The interplay between network dynamics and network 
topology is the key to understand collective behaviors on com- 
plex networks. Previously, many works have revealed how 
network topology affects the network synchronizability under 
the setting of coupled identical oscillators. For nonidentical 
oscillators, there are two factors affecting the development of 
GS on complex networks, i.e., network topology and hetero- 
geneity in the local dynamics. The effect of heterogeneity in 
the local dynamics, if they are chaotic, on the development of 
GS remains to be a challenging problem. In the two preceding 
sections, we have shown that for scale-free networks, GS typ- 
ically starts from the hubs, and then spreads to others nodes 
with relatively smaller degrees. In this case, heterogeneity in 
the degree distribution appears to be the main factor govern- 
ing the development of GS. An interesting issue is then that, 
for networks with homogeneous degree distributions, such as 
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FIG. 8: (Color online) Color-maps of d(e,id) and d(e,i r ) in the 
two-dimensional parameter space (e, i) for nonidentical Lorenz os- 
cillators on a regular network. The network has 100 nodes, and each 
node has 6 nearest neighboring connections. 
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FIG. 9: (Color online) Color-maps of d(e, id) and d(e, i r ) in the two- 
dimensional parameter space (e, i) for nonidentical Lorenz oscilla- 
tors on a small- world network. The network is obtained by rewiring 
a small part of connections on a regular network with 100 nodes 0]. 



regular, small-world, or certain modular networks, how does 
the heterogeneity in the local dynamics affect GS? In the fol- 
lowing, we study this question numerically. 

We first consider a regular network of N = 100 nodes. 
Each node in the network has kd = 6 connections to its nearest 
neighbors. As a special case, this is an exactly homogeneous 
network with a regular degree sequence. The local dynamics 
are those of Lorenz chaotic oscillators with different param- 
eter n randomly distributed in the interval [28, 38] for this 
example and for the following two examples. Figure 8 shows 
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FIG. 10: (Color online) Color-maps of d(s 1 i c i) and d(s,i r ) in the 
two-dimensional parameter space (e, i) for nonidentical Lorenz os- 
cillators on a clustered network. The network has 100 nodes which 
are evenly divided into 5 modules. Nodes within each module are 
fully connected. Any two nodes in different modules have probabil- 
ity p — 0.005 to connect each other. 



the color-maps of the local synchronization distances d(e,id) 
and d(e, i r ), where id and i r are the node indices arranged ac- 
cording to the decreasing node degree and decreasing values 
of parameter r^, respectively. For the Lorenz oscillator under 
the parameter setting in our study, we find that the larger the 
value of parameter r, the larger the largest Lyapunov exponent 
of the chaotic attractor. Thus, in Fig. 8, the index i r actually 
corresponds to the decreasing value of the largest Lyapunov 
exponent. Specifically, larger value of i r corresponds to a 
small value of T{ so that the corresponding local dynamics 
is less chaotic in the sense that its largest Lyapunov exponent 
has a relatively smaller value. Comparing Figs. 8(a) with 8(b), 
we find that, when transient behaviors are disregarded, nodes 
whose dynamics are less chaotic require smaller value of the 
coupling strength to be entrained with other nodes in the net- 
work. Entrainment of more chaotic nodes requires stronger 
coupling. Thus the development of GS on regular network is 
dominated by the local dynamics. 

We next consider the occurrence and development of GS 
on small- world networks. A representative example is shown 
in Fig. 9. We see from Fig. 9(a) that there is no apparent 
synchronization sequence of nodes according to the degree 
index i&. However, as can be seen from Fig. 9(b), GS starts 
from nodes with smaller values of the largest Lyapunov expo- 
nent. This demonstrates that for a small- world network, het- 
erogeneity in the local dynamics plays a dominant role in the 
development of GS, which is similar to the situation in regular 
networks. 

Lastly, we investigate GS in a type of modular network. 
A modular network is characterized by a number of sparsely 
connected subnetworks, each with dense internal connections. 
The modular structure has been identified in biological, social, 



and certain technological networks 112811 . For such a network, 
synchronization within each individual cluster can usually be 
achieved readily due to the dense internal connections, so the 
occurrence of global synchronization is particularly interest- 
ing Jl9ll . An example of the development of GS on modu- 
lar network is shown in Fig. 10. We see that global GS can 
be achieved despite the sparse inter-cluster connections. An 
interesting phenomenon is that, the development of GS does 
not seem to strongly depend on id or i r . Nevertheless, com- 
paratively, it still can be found that nodes with less chaotic 
local dynamics are easier to be entrained with smaller cou- 
pling strength. This can be explained by noting that the de- 
gree distribution is approximately homogeneous, as a result of 
the sparse inter-cluster connectivity. In addition, nodes with 
different values of r are uniformly distributed in different syn- 
chronous clusters. As a result, the effect of heterogeneity in 
the local dynamics is suppressed to some extent. 

Based on the above studies, we can conclude that the de- 
velopment of GS on complex networks is the result of interac- 
tion between network topology and local dynamics when they 
both have heterogeneity. For heterogenous networks, such 
as scale-free networks, the network topology plays a leading 
role; while for approximately homogeneous networks, such 
as small- world networks and modular networks, the local dy- 
namics is the dominant factor affecting the development and 
spreading of GS on networks. 



V. EFFECT OF DIFFERENT COUPLING STRATEGIES 
A. The linear coupling and the normalized coupling 

We have seen that the hub nodes behave as "seeds" to de- 
velop GS in scale-free networks. The phenomenon that syn- 
chronization starts from the hub nodes in a heterogeneous net- 
work is in fact quite general, which holds for both GS and CS. 
Qualitatively, this can be understood by regarding the inter- 
actions among nodes as a kind of "noisy" driving within the 
network. In particular, since the dynamics of each individual 
node is chaotic, the coupling term in Eq. (Q} for node i is 
also chaotic, which can be regarded effectively as correlated 
noisy driving. It has been known that correlated noise can in- 
duce and/or enhance synchronization [29], PS IH, and GS 
as well |4|]. That is, correlated noise of larger amplitude can 
cause oscillators under such driving to be better synchronized. 
For a heterogeneous network, hub nodes have significantly 
more links than other nodes. Since the drivings from differ- 
ent nodes are not entirely uncorrelated, the equivalent driving 
force to the hub nodes are generally stronger, facilitating their 
synchronization with other nodes in the network. 

However, in a recent publication studying GS on networks 
[27], it is reported that only for a kind of special scale-free net- 
work with tree-like structure, GS is observed to develop from 
the hubs and then gradually spread to other nodes in the net- 
work; while for usual scale-free networks, the heterogeneity 
of network seems to have little effect on the development of 
GS in the network. In attempting to find the reason that leads 
to different observations in the present work and Ref. |27|], we 



8 



notice that a different coupling strategy is used in Ref. l27h . 
Take a coupled one-dimensional map network as an example, 
the coupled system in Ref. [12711 is: 

x n+l = {1 _ £)/(a .n ) + L ^OijfW), (6) 

* 3 

where x™ is the state variable of node i at time step n, f(x) 
is the local map, and k{ is the degree of node i. This cou- 
pling scheme with each node in the network driven by the 
local mean field is essentially the same as the following nor- 
malized coupling strategy: 

Xi = F^x*) - |- y^q^-(xj - x^), (7) 

where notations are the same as that in Eq. (Q]) and Eq. ©. 
In this coupling scheme, the coupling strength for each node 
i is normalized by its degree. Therefore, the effective cou- 
pling strength each node received from its network neighbors 
is of the same order of magnitude, regardless it is a hub or a 
node with very small degree. Obviously, due to the normaliza- 
tion, the effect of network topology on the development of GS 
has been almost suppressed. As a consequence, there should 
be no distinct difference among the coupling strength when 
nodes in the network achieve GS. Quantitatively, we can find 
that in order to obtain the same effective coupling strength, 
system © should require (k) times larger coupling strength 
than system (Q]). Here, (k) denotes the average degree of the 
network. Thus on an average, it can be expected that system 
© needs (k) times larger coupling strength to achieve GS in 
network compared with system (Q]). 

The above heuristical considerations can be further analyti- 
cally explained. By examining the coupled system (Q]) and its 
corresponding auxiliary system (|2]). Letting Ax^ = x • — x^ 
and subtracting Eq. (Q]) from Eq. 0, we obtain 

Axi = Fi(xJ) -Pi(xi) - e^OijAxi (8) 

3 

« [DF(x^) — ekil] • Axi, 

where I is the M x M identity matrix with M the dimen- 
sion of local dynamics. We see that, for node i, the effec- 
tive coupling strength is proportional to its degree k\. For 
a fixed value of e, the coupling between a hub node and its 
counterpart in the auxiliary system can be significantly larger 
than the coupling for nodes with smaller degree, leading to 
"earlier" synchronization between the hub nodes in the orig- 
inal and in the auxiliary systems. Therefore, for usual linear 
coupled oscillator system, the general observation is that, in a 
complex network with heterogeneous degree distribution, the 
set of hub nodes provides a skeleton around which synchro- 
nization is developed. The above analysis is also suitable for 
coupled system ©. In this case, there will be no factor k{ as 
in Eq. ([5]). This implies that the effect of network topology on 
the development of GS in networks is eliminated. 

In the following of this section, we further present two ex- 
amples where the oscillators in a scale-free network are inter- 
acted with the normalized coupling strategy. The first exam- 
ple is coupled Lorenz oscillators described by Eq. ©, and the 
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FIG. 11: (Color online) Color-maps of d(s,id) in the two- 
dimensional parameter space (e,i), characterizing the development 
of GS for 300 identical Lorenz chaotic oscillators (a), and 300 iden- 
tical Logistic maps (b) on a scale-free network. The network is the 
same as in Fig. 1. Compared with Fig. 1, the effect caused by the 
normalized coupling strategy is obvious. 



second example is the following coupled Logistic maps: 

*r +1 = / w) - f - /(*?)]• w 

3 

For both cases, the local dynamics are identical, i.e., r^=28 
for all Lorenz oscillators and f(x) = 4x(l — x) for all Lo- 
gistic maps. Figure 1 1 illustrates the development of GS for 
the above two systems on a scale-free network. The effect of 
different coupling strategies on the development of GS can be 
verified by comparing Fig. 11 with Fig. 1. For system (Q]), 
GS first occurs on the hubs and then gradually spread to other 
nodes; while for systems © and ©, GS almost simultane- 
ously takes place on the hubs and the other nodes. Actually, 
there is no significant difference among all nodes on networks. 
By examining Fig. 1 and Fig. 1 1(a), it is found that the giant 
GS cluster forms at e « 0.5 and e w 2, respectively, in the 
two cases. In particular, the latter is approximately four times 
larger than the former. Note that the scale-free networks we 
used in Fig. 1 and Fig. 1 1(a) are BA model with mo = m = 4, 
so the average node degree (k)=4. This is consistent with the 
theoretical analysis above. 

B. Coupled hybrid oscillators: competition of local dynamics 

In Sec. IV, we have studied how network topology com- 
petes with local dynamics to dominate the development of GS 
in complex networks when there are heterogeneity in both net- 
work topology and local dynamics. In fact, for coupled non- 
identical oscillators on complex networks, competition may 
also occur among local dynamics if the parameter mismatches 
of local dynamics are large enough, or the local dynamics are 
physically different. By adopting the normalized coupling 
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FIG. 12: (Color online) Color-map of d(s,id) in the two- 
dimensional parameter space (e,i), characterizing 300 hybrid 
Lorenz chaotic oscillators on a scale-free network. The network is 
the BA model with ra = m = 3 Q. 




FIG. 13: (Color online) Color-map of d(e,id) in the two- 
dimensional parameter space (e, i), characterizing 300 hybrid 
Lorenz and Rossler oscillators on a scale-free network. The network 
is the BA model with mo = m = 3 fl3h . 



strategy, we can conveniently investigate this issue. In the 
following, we present two examples showing how GS is de- 
veloped for coupled hybrid oscillator system in networks. By 
hybrid we mean that the oscillators in the network can be clas- 
sified into different types which are parametrically different or 
physically different. In the first example, we consider a sys- 
tem coupled with hybrid Lorenz oscillators. To be specific, 
we arbitrarily select 5% Lorenz oscillators in the network and 
make their parameter n be randomly distributed in the inter- 
val [30, 40]. For the rest oscillators, they are the same with 
r^=28, which are significantly smaller than that of the 5% os- 
cillators. In Fig. 12, the development of GS for such a hybrid 
system is illustrated on a scale-free network. It is seen that 
most oscillators in the network are entrained at e w 2. How- 
ever, a small number of oscillators require significantly larger 
coupling strength to be entrained. A careful examination of 
the locations of these nodes reveals that they just correspond 
to the 5% oscillators with larger r values. In the second ex- 
ample, the coupled hybrid system consists of two kinds of os- 
cillators which have different dynamical equations. Similar to 
the first example, 5% nodes are randomly selected to be the 
identical Rossler oscillators, while the rest 95% percent nodes 
are the identical Lorenz oscillators with r^=28. The dynamical 



equations of the Rossler oscillator read: 

P(xi) = [-( yi + 2a), x % + 0.2^, 0.2 + zi{xi - 5.7)] T , (10) 

where x$ = (x^ y%,Zi) are the state variables of the Rossler 
oscillator. In Fig. 13, we plot the local distance of GS vs the 
coupling strength for all nodes in a scale-free network. It is 
clearly shown GS can be achieved for such a coupled hybrid 
system. The development of GS in network has two distinct 
stages: the Rossler oscillators are much easier to be entrained 
compared with the Lorenz oscillators. The latter needs quite 
larger coupling strength to achieve GS. Note that the largest 
Lyapunov exponent of the Lorenz oscillator is proportional to 
its r value, and the Rossler oscillator has much smaller Lya- 
punov exponent than that of the Lorenz oscillator, we can con- 
clude from the above two examples that GS usually develops 
from the nodes in the network where the local dynamics are 
less chaotic in the sense that the largest Lyapunov exponents 
of the local dynamics have relatively smaller values. Although 
our numerical simulations are carried out on a scale-free net- 
work, we believe this conclusion can still hold for other net- 
work topologies when the normalized coupling strategy is ap- 
plied. 

VI. SUMMARY 

In summary, we have investigated the occurrence and de- 
velopment of GS on various complex networks, including 
scale-free networks, small- world networks, random networks 
and modular networks. It is shown that GS generally takes 
place in such networks for both coupled identical oscillators 
and nonidentical oscillators. For coupled identical oscillators, 
we reveal a general bifurcation path to global CS, i.e., non- 
synchronization PGS -> GS -> PCS -> CS. Prior to the on- 
set of global GS (CS), there exists a coupling regime of partial 
GS (CS) characterized by the occurrence of a cluster of syn- 
chronous nodes. In the transient regimes, the size of the syn- 
chronous cluster usually oscillates with the increase of cou- 
pling strength, showing a dynamic scenario when the coupled 
system is approaching global synchronization. We find that 
the development of GS on complex networks depends on both 
network topology and local dynamics. Moreover, the specific 
coupling strategy also plays an important role during the evo- 
lution of synchronization. We show that under the usual linear 
dissipative coupling scheme, for heterogeneous networks, GS 
generally starts from a small number of hub nodes, and then 
spread to the rest nodes in the network; while for homoge- 
neous networks, GS usually starts from the nodes whose local 
dynamics are less chaotic in the sense that the largest Lya- 
punov exponents have relatively smaller values. We further 
show that the effect of network topology on the development 
of GS can be suppressed if the coupling strengths of nodes 
in the network are normalized. Under such coupling scheme, 
the development of GS is essentially determined by the com- 
petition among different local dynamics. We demonstrate that 
GS can also occur in coupled systems with hybrid oscillators 
on complex networks, and the development of GS has distinct 
stages corresponding to physically different local dynamics. 
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We note that, while there are extensive works on synchro- 
nization in complex networks, prior to this work the bifur- 
cation path toward global CS via GS, the dynamical features 
of the development of GS, especially the interplay between 
the network topology and local dynamics on typical complex 
networks have not been investigated. Our work reveals that 
on complex networks, coupled oscillators may present fun- 
damentally different synchronization regimes which deserves 



further study. 
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